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Abstract 

We calculate the leading nonperturbative contributions to the effective Wil- 
son coefficient of the four-fermion operators arising from the QCD penguin, 
and we demonstrate how the usual perturbative one loop contribution is aug- 
mented by nonperturbative condensates. These corrections, which are ob- 
tained by quark and gluon condensate insertions into the quark loops, enter at 
the next-to-leading logarithm precision. Our results indicate for the charmed 
quark loop that the gluon condensate contribution is quite sensitive to the 
momentum transfer to the quark- antiquark pair. 
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Nonleptonic B decays to light hadrons have recently received a lot of theoretical attention. 
This is mostly due to the measurements by the CLEO Collaboration of a number of B — > 
hifi2 decays, where h\ and h 2 are light hadrons such as n, K, rj' and uj []!]-§]. For example, the 
larger than expected experimental results for the branching ratio of the inclusive B —>■ X S 1]' 
and its exclusive counterpart B — > Krj' have led to various theoretical explanations both 
within and beyond the Standard Model On the other hand, these decays could provide 

the first evidence for CP violation in B decays in the near future ||. 

The standard theoretical framework to study the nonleptonic B decays is based on the 
effective Hamiltonian approach. In this way, the heavy degrees of freedom, i.e. top quark 
and W ± gauge bosons, are integrated out and the matrix element can be written as 

< hMH C K\B >= ^E^CKMC,(/i) < hMOi{n)\B > (1) 
v '1 i 

where {0{\ is the complete set of operators relevant to B — > hih 2 decay and Ci(fi) are 
the corresponding Wilson coefficients evaluated at the renormalization scale fi [K|. These 
coefficients have been calculated up to the next-to-leading logarithm (NLL) precision in the 



strong coupling constant a s IjlUJEq] . At this precision, the matrix elements on the RHS of 
that equation should be evaluated at the one loop level in order to have a renormalization- 
scheme independent result for the LHS of the eq. (1) |TTi|l3 . 



In this work, we concentrate on the four-Fermi operator subset of {Oi}, i.e. i = 1..6, 
which is relevant for the hadronic B decays. The first two operators 0\ and 2 are the 
current-current opertors: 

Oi = qofi^LUaUp^^Lbp , 2 = q a ^Lupu^^Lb a , (2) 

and the other four are the so-called QCD penguin operators: 

O3 = q a YLb a q'pipLqp , A = q a YLb p ^ q^^'a » 

q> q> 

C 5 = QaYLb a V'l3ln R( l'l3 , °e = qal^Lbp $i7/^a • ( 3 ) 

q> q> 

where q = s or q = d for b s or b — > d transitions, respectively, and where L and R are 
the left and right handed projection operators, i.e. L = 1 — 75 and R — 1 + 75. Two more 



current-current operators which are obtained by replacing u quark by c quark in eq. (2) 
should also be added to the above list. 

In NLL precision, the matrix element of the operators is to be evaluated at one loop 
level, as mentioned above. These one loop matrix elements in turn can be written in terms 
of tree level matrix elements with loop effects included in the corresponding effective Wil- 
son coefficients C? s ||TTlJl3|| . There are two types of relevant one loop diagrams, the vertex 



corrections, where a gluon connects two of the outgoing quark lines, and penguin diagrams 
where a quark-antiquark pair forms a loop and emits a gluon which in turn produces a 
quark- ant iquark pair. Our focus in this paper is on the latter type of corrections (Fig. 
1) since they could have a potentially significant effect on branching ratios and CP asym- 
metry in two body hadronic B decays |13|]l4}| . In fact, the absorptive portion of penguin 



diagrams can provide the strong phase necessary for direct CP violation. This one loop cor- 
rection is analogous to the so-called continuum contribution in the dileptonic rare B decays 
B — > X s {X^)t^ 'i~ ', where one has to deal with a lepton-antilepton pair, as opposed to a 
quark-antiquark pair \IE\. For the dilepton case, we also have significant nonperturbative 



contributions from intermediate vector mesons which are quite sensitive to the momentum- 
transfer to the dileptons. Motivated by the above analogy, we embark on supplementing 
Cf with a corresponding set of nonperturbative contributions. 

In this work, we calculate the leading nonperturbative corrections to the perturbative 
evaluation of the penguin diagrams. Ignoring the vertex corrections, only the QCD penguin 
operators receive corrections from one loop diagrams of Fig. 1 |13j : 

cf = c 9 -\%(c t + c p ) + ... , 
cf = c A + 1 -^(c t + c p )+ ... , 

cf = C 5 — (C t + C p ) + ... , 

Cf = C 6 + 1 -^(C t + C p )+ ... . (4) 

C t and C p are the loop corrections corresponding to the current- current (tree) and penguin 
type four-fermion operator insertion in Figure 1, respectively. Using naive dimensional 



regularization(NDR) along with the MS renormalization scheme, one obtains 
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where the function AFi(z) is defined as 

AFAz) = -4 / - a;)/n [1 - zx(l - x)\ 

Jo 



(6) 



We note that the effective coefficients in NLL depend on the momentum-transfer k 2 as well 
as the internal quark mass m q . In fact, for z > 4, AFi(z) develops an imaginary part 
which has been taken as the strong phase necessary for generating direct CP asymmetry 
|T3| . Note that the light internal quarks produce the dominant contributions; consequently, 
the larger factor C t (C\ is much larger then Cz...q) is more significant for b — > d transitions 
in which the ratio of CKM factors for tree and penguin is not too small. Nevertheless, 
one can not have a complete picture of the NLL corrections without taking into account 
nonperturbative effects which, as in the case of the intermediate resonance contribution to 
B — > X s {Xd)C + d~ i can prove to be significant. To investigate this possibility, we calculate 



the leading nonperturbative contributions to C, 



3.. 6 by inserting quark and gluon 



condensates, which are believed to be responsible for the chiral-symmetry breaking and 
confinement properties, respectively, of the QCD vacuum, into one-loop penguin diagrams 
(Figs. 2,3). 

QCD condensates characterize the nonperturbative content of the QCD vacuum, and 
therefore, condensate insertions should give us a hint of the leading nonperturbative con- 
tributions to the effective Wilson coefficients. To demonstrate our method, we present a 
detailed calculation of the quark condensate result. For the gluon condensate contribu- 



tion, the necessary expressions are extracted from previously determined |16j condensate 
contributions to the vector correlation function. 



To obtain the quark condensate contribution, we proceed by replacing the usual pertur- 
bative internal quark propagator in Fig. 1 with the full quark propagator S(p): 

S(p) = S p (p) + S NP (p). (7) 

The nonperturbative contributions to the quark propagator can then be written as 

iS NP {p) = J d 4 xe ip - x < fi| : (0) : |0 > , (8) 

where \Q > is the nonperturbative QCD vacuum. The nonlocal vacuum expectation value 



(vev) in eq. (8) can be expanded in terms of local condensates fl7| . To ascertain how 
the leading nonperturbative QCD contributions affect the coefficients C[ s , only the lowest 
dimensional (quark and gluon condensates) components of the expansion are considered 
below. However, we note that the nonlocal vev of two quark fields does not contain a gluon 
condensate component. The quark condensate projection of the nonperturbative quark 
propagator is taken from Ref. ||18| : 



iS NP (p) = (2n)\p' + m q )F(p) (9) 
where the Fourier transform of F(p) is expressed in terms of Bessel function 

/^Fo.)=- y Ji(, ^f ' . no) 

with the following important property: 

(p 2 - m 2 q )F(p) = . (11) 

The dimension-3 quark condensate < qq > is the vev of the normal ordered local product 
of quark and antiquark fields, i.e. 

< qq >=< Q\ : f (0)^(0) : \Q > . (12) 

Using the Feynman rule of eq. (9), one can obtain the nonperturbative < qq > contribution 
to the effective Wilson coefficients in a straightforward manner. The relevant Feynman 



diagram is illustrated in Fig. 2, where the nonperturbative quark propagator 5' NP is depicted 
by a disconnected line with two dots. Here we concentrate on the loop portion of Fig. 2 
which differs from the perturbative case of Fig. 1. Aside from the color factor, the vector 
current correlation function can be written as 



,i Tr[{$- # + m)7 M (2$ + m)7„]F(p) 

! h ! Z I ( \ 

'" 1 (p — k) — m? + ie 



n<*>(A0 = 2 d 4 P 



(13) 



where the factor 2 in front is due to two possible insertions of the nonperturbative quark 
propagator in the fermion loop. By contracting p^p v into n<J 9> and using the identity 

d 4 pk.pF(p) = i\im-^- J d 4 pe- iik - p F(p) 

= l <( ii > lim j jiK^j = o 



_ (14) 
6m q r-A)dC m^Vk 2 

where the second line is derived via eq. (10), one can show explicitly the transversality of the 
correlation function in eq. (13). Consequently, one finds upon contracting g^ v into n<J 9> 
and imposing the on-shell constraint (11) that 

F(p) 



U <?q>ll Jk) = 24(2 



m~ + k z 



d*p-. 



k 2 — 2k.p + ie 

The integrals appearing in eq. (15) are evaluated as follows [15]: 

< qq > 
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(16) 



where the final line is derived by utilizing a tabulated integral lj| . Using the transversality 
of II<J 9> , the result for the RHS of (13) is easily obtained: 

< qq > 



3m^ 



1 



1 



2m 2 q 



1 



Am 2 
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(17) 



k 2 J V k 2 

This result is, in fact the < qq >-contribution to the vector-current correlation function 



derived in ref. [IB]. The corresponding gluon condensate contribution to the vector-current 
correlator, which may be derived by several different methods [0, is given by [|l(| 
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where 



z = k 2 /m 2 



(19) 



and where |T6| , pO 
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It is evident from (18) and (20) that this contribution diverges sharply at z = 4 pT |. 

Consequently, we replace eq. (6) with the following < qq >- and < G 2 >-augmented def- 
inition for A.Fi in order to include the leading nonperturbative quark and gluon condensate 
contributions to the effective Wilson coefficients (eq. (4)): 



-4 / dxxil — x)ln [1 — 2x(l — x) 
Jo 
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(21) 



where the second term is written in terms of the renormalization group (RG) invariant 
combination < m q qq >, and where the factor a s in the third term is absorbed in the gluon 
condensate order parameter to achieve approximate RG-invariance. 

In the decays B — > h\h,2, where hi and fi2 are light mesons, the characteristic range of 
k 2 is argued to be p2[ 



2 2 
m b < , 2 < m b 
~ K ~ . 

4 2 



(22) 



For light quarks (z = k 2 /m 2 q » lj, the nonperturbative terms in eq. (21) can be approxi- 
mated by 



^ ^nonperturbative 



lQn 2 < m q qq > 2ty < a s G 2 



> 



(23) 




Inserting the phenomenological values < qq >~ (— 0.25) 3 GeV 3 and < aG 2 >~ 0.045 
GeV 4 , we observe that the leading nonperturbative corrections are around three orders of 
magnitude smaller than the perturbative contribution, and can therefore be safely neglected. 

For the heavier charmed quark, however, the situation is more interesting because of the 
presence of a singularity at z = 4 in the gluon-condensate term: the range AGeV 2 ~ k 2 ~ 
lOGeV 2 suggested by (22) is quite likely inclusive of z — 4, (i.e., of k 2 = Am 2 ) given the 
present 1.1-1.4 GeV range of values for m c ||23|| . Consequently, there is a genuine possibility 
that the gluon-condensate contribution to the penguin amplitude could be comparable to the 
purely perturbative correction, depending on which value of k 2 one chooses. This sensitivity 
to the momentum transfer k 2 is also important for investigating direct CP asymmetry in 
hadronic B decays, as the gluon-condensate contribution to (21) develops an imaginary part 
which dominates the strong phase of the effective Wilson coefficients as z — > 4 from above. 

For a physical interpretation of the above results, we again use the analogy with the 
dileptonic B decays B — > X s £ + £~ . These decay modes receive a perturbative charmed-quark 
loop correction (continuum contribution) as well as nonperturbative contributions when cc 
forms intermediate resonances such as J/ip, if)',--- ||15|| . The latter can be incorporated by 
modelling a Breit-Wigner form for the resonance propagator in the long distance amplitude 



where k is the momentum transfer to dileptons, and where M and T are the mass and total 
decay width of the intermediate resonance. The resulting dileptonic differential decay rate 
has peaks at k 2 ~ M 2 , and the nonperturbative resonance contribution indeed dominates 
the branching ratio of this mode. 



PI 




1VC 



1 



(24) 
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One can make a similar charmonium-resonance interpretation for the singularity at k 2 = 
Am 2 in the < G 2 > contribution to the effective Wilson coefficients of the condensate- 
augmented penguin amplitude, particularily if we identify the k 2 > 4m 2 kinematic threshold 
for the imaginary part with the onset of a "physical" intermediate state. It must be noted, 
however, that such a resonance couples directly to a single gluon (Fig. 3), necessarily 
implying that the z = 4 singularity, if indeed a resonance, must be associated with the 
leading contribution of a weakly bound colour-octet cc state [pq] . 

In conclusion, we have calculated the nonperturbative quark and gluon condensate con- 
tributions to the effective Wilson coefficients of the QCD penguin four-fermion operators, 
supplementing the one loop (continuum) qq correction which appears at NLL precision with 
its nonperturbative counterparts. We have pointed out how incorporation of the gluon con- 
densate in the charmed quark loop results in a contribution which is highly sensitive to 
momentum transfer k 2 . Since k 2 is not in-and-of itself a physical observable in two body 
hadronic B decays, one cannot exclude the possibility of significant nonperturbative contri- 
butions to both the magnitude and the strong phase of the effective Wilson coefficients of 
the four-Fermi operators O^...Oq. 
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Figure Captions 



Figure 1: Perturbative one loop correction to the effective coefficient of the four-fermion 
operators. 

Figure 2: Nonperturbative quark condensate contribution to the effective coefficient of the 
four-fermion operators at the one-loop level. 

Figure 3: Nonperturbative gluon condensate contribution to the effective coefficient of the 
four-fermion operators at the one-loop level. 
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